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ABSTRACT 



In this thesis, a new scheme of encoding orthogonal MFSK signals with convolutional codes 
is proposed. Instead of separating the code and the signals as is done conventionally, here the 
m outputs of the convolutional encoder of rate k/m are used to select the M = 2 m orthogonal 
FSK signals. This combined coded system is called trellis-coded modulation MFSK 
(TCM/MFSK). 

The free distance of TCM/MFSK codes is found computationally. Asymptotic coding gains 
(ACG) of several TCM/MFSK systems relative to uncoded MFSK systems having the same 
bandwidth efficiency (bits/ s/Hz) are calculated. 

Both analytic descriptions and natural mapping implementations of TCM/4-FSK with rate 1/2 
and TCM/8-FSK with rate 2/3 for several constraint lengths of the convolutional codes are 
given. The analytic description of TCM/MFSK is also obtained computationally. 
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I. INTRODUCTION 



The model of the digital communication system under consideration is shown in 
Figure 1.1. The input to the convolutional encoder, {aj, aj k = 0, 1 is a sequence of k 
independent bits. The channel encoding operation is combined with M-ary Frequency 
Shift Keying (MFSK) modulation as in Trellis Coded Modulation (TCM) schemes. A 
convolutional encoder is used to generate the underlying code. MFSK orthogonal 
signaling is used to achieve an acceptable performance with a minimum power 
requirement. 

The set of orthogonal signals used is denoted as: 

{S x ( t)}, *=1,2, M (1.1) 

These signals are generated by a frequency synthesizer which is driven by a set of real 
numbers 



n<t) 




Figure 1 . 1 System Model 



1 



{x}, x=l,2, M 



( 1 . 2 ) 



These real numbers have a specific relationship with the incoming bit stream {a,} and the 
memory bits of the encoder. For an encoder that has 2 V states, the constraint length of 
the convolutional code is defined as v. The constraint length v is equal to the memory 
of the encoder (shift register stages). 

The M = 2 m distinct symbols are represented by M orthogonal sinusoidal 
waveforms as: 

S x {t) = A cos [2 ir(f c + (x - 1) A fit ] , 0 <t<T (13) 

where 

x = 1, 2, M 

T: signaling interval 
f e : carrier frequency 
Af: the width of the frequency slot. 

The orthogonal signal set is characterized by equal signal energy E given by: 

= (1-4) 

The orthogonality condition is satisfied when the cross-correlation coefficient obeys the 
relationship 



2 



T 






0 ifi*j 

E if H 



(1.5) 



MFSK is generated by subdividing a frequency interval into M distinct frequency slots. 
Each slot has a width Af. The minimum frequency separation between adjacent signals 
is obtained from (1.5) [Ref. 1] 



Af = 

J mm 



T 



( 1 . 6 ) 



Thus, the minimum one-sided bandwidth occupied by these M orthogonal signals is: 




(1.7) 



It is assumed that the transmission medium introduces zero mean, additive white 
N 

Gaussian noise (AWGN) n(t) with power spectral density . The received signal is 

2 

corrupted by AWGN as: 

r x (t)=s x (t)+n{t) (1.8) 

This signal, r x (t), is demodulated by a bank of either matched filters or correlators. 
When the initial phase can be estimated by the receiver, the demodulation is coherent, 
otherwise it is noncoherent. 

The analog outputs (u x (t)}, x = 1, 2, ... M of the demodulator are then sent to the 
Viterbi decoder which performs a soft decision decoding of the MFSK signal. The 
sequence {dj} at the output of the Viterbi decoder is the maximum-likelihood estimate of 
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the input sequence {a,}. In this thesis, only coherent TCM/MFSK systems are 
considered. 

The remaining chapters in this theisis are organized as follows. In Chapter n, an 
overview of convolutional coding and TCM/MFSK are presented. The results from the 
computer program written to find the free distance of TCM/MFSK are given. Asymptotic 
coding gain (ACG) is defined and a list of ACGs is also given. Chapter HI deals with the 
design of the encoder for TCM/MFSK. Both the analytic expressions and the natural 
mapping rules for certain schemes are examined. Although all the signals in the M-ary 
orthogonal constellation have equal distances from each other, it is shown that the 
transmitter complexity depends on the signal assignment. Basic rules are given for the 
minimal transmitter complexity. In Chapter IV, several encoders are designed on a step 
by step basis. TCM/MFSK is further discussed and some conclusions are drawn in 
Chapter V. 

The computer program to find the free euclidean distance of TCM/MFSK is given 
in Appendix A. Another program to find the analytic expression of TCM/MFSK is given 
in Appendix B. 
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II. OVERALL REVIEW 



In 1948, Shannon showed that the main limitation to transmission over a noisy 
channel is set by a quantity called channel capacity. Capacity of a noisy channel is the 
largest rate at which information can be transmitted reliably. In other words, if the data 
source rate is less than the channel capacity, proper encoding and decoding techniques 
enable us to communicate over a noisy channel with any arbitrary error rate. Otherwise, 
reliable communication is not possible. 

Codes are divided into two categories according to their function. Codes for source 
encoding attempt to reduce the data source rate at the price of tolerable distortion. Codes 
for channel encoding turn a noisy, constrained channel into a reliable, unconstrained 
channel. The application of channel encoding is called Forward Error Correcting Coding 
(FEC). Most communication systems use FEC coding systems. FEC coding is generally 
achieved by redundancy adding and noise averaging. Codes are classified according to 
their structure as block codes and convolutional codes. 

In this thesis, channel encoding and MFSK modulation operations are combined 
together to take advantage of soft decision decoding. Convolutional coding is used for 
this purpose. It may be possible to use another code as the outer code to improve 
performance further. 

The basic distinction between convolutional codes and block codes is that 
convolutional codes do not have a fixed block length. These codes are generated by 
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passing the input sequence through a system of linear shift registers and modulo-2 
adders. A convolutional code is characterized by a ratio called the code rate. The code 
rate is the ratio of the number of input lines to the number of output lines of the encoder. 
Thus, the code rate is 




m 



where k is the number of inputs and m is the number of outputs of the encoder. The 
performance of the encoder depends strongly on the number of states 2 V where the 
constraint length v is equal to the memory of the encoder. 

The basic utility of convolutional codes comes from the fact that these codes can 
be decoded by the Viterbi decoding algorithm. The Viterbi algorithm is an efficient way 
of sequence demodulating in the presence of noise [Ref. 2]. The Viterbi decoding 
algorithm is equivalent to maximum-likelihood decoding and thus optimum for equally 
likely messages [Ref. 3], 

When the encoding and modulation are done separately, the basic parameter 
affecting the code performance is the Hamming distance, which is the number of places 
in which any two sequences differ. If the encoding and the modulation are combined into 
one step, the Euclidean distance governs the performance of the coded scheme. The 
Euclidean distance between two sequences c and d is defined as 

^(c,</)=[£m \c r d t \ 2 ] m ( 2 - 2 ) 

where c, and d t are the elements of c and d, respectively. 
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At the receiver end, if the Viterbi decoder deals with the quantized demodulator 
outputs which are the best estimates of the transmitted signal (i.e., the number of 
quantization levels is 2), it is called hard-decision decoding. If the Viterbi decoder works 
on the sequence of information on which a final decision can be made about the 
transmitted signals (i.e., the number of quantization levels is more than 2), the method 
is called soft-decision decoding. Soft-decision decoding has more coding gain than 
hard-decision decoding [Ref. 1]. 

Trellis coded modulation (TCM) proposed by Ungerboeck combines encoding and 
modulation to design a code for a particular geometry of the signal constellation. Thus, 
the Euclidean distance becomes the main parameter and soft-decision decoding is always 
used. The trellis encoder is composed of two parts. The first part is the convolutional 
encoder with v bits of memory. This determines the number of the states 2 V of the 
encoder. The second part is the mapper. This part is memoryless. Coded digits select the 
signal from the constellation and move the encoder to the next state. Uncoded digits only 
select signals causing parallel transitions. The mapper maps the incoming coded and 
uncoded digits to the M signals. The mapping rule is defined by a function which must 
be nonlinear in order to achieve a coding gain [Ref. 4]. 

The asymptotic coding gain is the ratio of Eb/N 0 required for coded operation 
relative to that required for uncoded operation with the same bandwidth efficiency in the 
limit of large E^/N,,. Although actual coding gain is a more realistic indicator of system 
performance, it is significantly more difficult to obtain than is ACG, and ACG does give 
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an indication of the performance that can be expected for AWGN channels. The ACG 
of TCM/MFSK is defined as: 

ACG = lOlo z(Xd\ e J2E) for TCM/4-FSK (2.3a) 

ACG = 101og(l4 £ /2 E) for TCM/S-FSK (2.3b) 

where 

d free : free Euclidean distance of the coded scheme (the smallest of the Euclidean 
distances between any two coded paths) 

E : average energy of the coded signal 

Decoding is done using the Viterbi algorithm. It is basically a search of the most 
likely path through the trellis. Sometimes the chosen path may not be the transmitted path 
due to the effect of noise. Then an error event is said to have taken place. An error event 
is the part of the path which diverges from the correct path and remerges at a later time. 
The free distance of a TCM scheme is the minimum Euclidean distance of all error 
events. Thus, by definition, free distance plays an important role in determining the 
performance of any TCM scheme. 

A general analytic solution for the free distance can not be obtained; hence, it is 
found computationally. A computer program is given in Appendix A to find the free 
distance of any TCM/MFSK. When the signal set is orthogonal, the distance between any 
pair of M signals is: 
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dy = \JlE for all ij 



(2.4) 



This is also the minimum distance of any orthogonal signal constellation. In the computer 
program this distance is normalized to unity regardless of the number of the signals in 
the constellation. Thus, 



The number of states in the trellis diagram is 2 V . The number of the transitions that 
diverge from a common state or remerge into a same state is 2 k . Since the code rate is 
r = k/m with rn_>_k + l,M = 2 m signals are assigned to 2 k+v transitions. The figure 
of merit to compare different coding schemes is ACG given by (2.3). The reference 
uncoded constellation is the one which has the same bandwidth efficiency as the coded 
scheme. Convolutional coding does not effect the spectrum of the transmitted signal 
[Ref. 5]. Hence, the bandwidth efficiency R/B may be taken as a reference for 
comparison purposes. Here R is the information rate in bits per second and B is the 
minimum one-sided bandwidth in Hz given by (1.7). For an uncoded M-ary orthogonal 
signaling, where M = 2 k ’ the information rate is 



dy = yJlE = 1 for all if 



(2.5) 




( 2 . 6 ) 



where k’ = log 2 M and T is the symbol duration. Also B = M/T, hence: 
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(2.7) 



R = = 

B M 2* 

When a convolutional encoder with rate r = k/m is used, the information rate is R = 
k/T, but the bandwidth is expanded to B = 2 m /T. Thus, the coded bandwidth efficiency 
is 

* = A (2.8) 

B 2 m 

ACG is calculated for the coded schemes for which (2.7) and (2.8) are equal. For 
example, a rate Vi TCM/4-FSK with k = 1 , m = 2 has a bandwidth efficiency of V* 
b/s/Hz. Uncoded 16-FSK with k’ = 4 has the same bandwidth efficiency. 

A search with the program given in Appendix A shows that the squared free 
distance of the best codes is 



4*“ - i + 1 (2.9) 

2E 



where L = |_ v /kJ , and L* J denotes the integer part. The results are tabulated in Table 
2.1. Note that the bandwidth efficiency for all TCM/4-FSK and TCM/8-FSK schemes 
is !4 b/s/Hz. 

The program input consists of a 2 V x 2 k matrix N. The rows of N are the transitions 
from the states of the trellis respectively. An example of a trellis diagram and the 

corresponding matrix N are given in Figure 2.1. The states are labeled as 1, 2, , 2 V . 

The first row of N consists of the labels of the destinations of the transitions from the 
first state. The second row for the transitions from the second state, so on. 
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1 



1 



2 



3 



4 



5 



6 



7 



8 




2 

3 

4 

5 

6 

7 

8 



N= 



12 3 4 

5 6 7 8 
12 3 4 
5 6 7 8 
12 3 4 
5 6 7 8 

12 3 4 

5 6 7 8 



Figure 2.1 A Trellis Diagram and its Matrix 
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TABLE 2.1 ASYMPTOTIC CODING GAINS OF TCM/4-FSK AND TCM/8-FSK 
WITH RESPECT TO UNCODED 16-FSK 



Code Rate 
(r) 


Modulation 

(MFSK) 


Constraint 
Length (v) 


2E 


ACG (dB) ! 


1/2 


4-ary 


2 


3 


-1.25 


1/2 


4-ary 


3 


4 


0.0 


1/2 


4-ary 


4 


5 


0.97 


1/2 


4-ary 


5 


4 


1.76 


1/2 


4-ary 


6 


7 


2.43 


1/2 


4-ary 


4 


4 


3.01 


1/2 


4-ary 


8 


4 


3.52 


1/2 


4-ary 


4 


10 


3.98 


1/2 


4-ary 


10 


If 


4.39 


2/3 


8-ary 


2 


2 


0.0 


2/3 


8-ary 


3 


2 


0.0 


2/3 


8-ary 


4 


3 


1.76 


2/3 


8-ary 


5 


3 


1.76 


2/3 


8-ary 


4 


4 


3.01 


2/3 


8-ary 


7 


4 


3.01 


2/3 


8-ary 


8 


5 


3.98 


2/3 


8 -ary 


9 


5 


3.98 


2/3 


8-ary 


10 


6 


4.77 
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m. ENCODER DESIGN 



The function of the encoder is to map the input binary stream into M orthogonal 
signals. This may be viewed as a sliding window method of encoding a binary data 
stream into a sequence of M real numbers. If the number of input bits into the encoder 
is k and the code constraint length is v, then the channel input S x (t) will depend on n = 
k + v bits. The v bits in memory determine one of 2 V states of the encoder while 2 k 
transitions are associated with k input bits. 

The MFSK waveform is given by (1.3) and it is the output of a voltage controlled 
oscillator (VCO) driven by a set of real numbers (voltage levels). Hence, the encoding 
operation may be implemented by mapping n bits into the set of real numbers {x}, x=l, 
2, ... M. Because x is a real-valued function of n bits, it may be written as a sum of 
products of input bits ^ where % E{0, 1} or b ; where b ( €{-1,1} and b; = l-2aj 
[Ref. 6]. Thus, 



n 



n 



*Q> v b v b n ) = d o +Y, d t b t + dybfi. 

i = 1 l V~ 

J>i 

n 

+ i^fi = i dijhbfijbb + d l2n b l b 2 b n 

h >j > i 



(3.1) 



where 



x (b„ b 2 , bj E {x}, x = 1, 2, M and 4, d^, d, jh , 
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are the coefficients that we want to find, and b„ b p b,, €{-1,1}. 
Following the discussion in [Ref. 4], we note that there are 





(3.2) 



coefficients. Let x denote a column vector with 2" elements which are all the values that 
x can assume. Also, let d denote a column vector whose elements are the unknown 
coefficients. Let B be a 2 n x 2“ matrix where each row represents the 2“ values taken by 
all the products of b;’s in (3.1). Thus, (3.1) may be written in matrix form as follows: 

x=Bd (3-3) 



where 



*(U,1, 1) 

*(- 1 , 1 , 1 , 1 ) 



x = 



*(- 1 ,- 1, -1 




2" xl 



and 



B k = [1 b, b 2 b 12 b 12 J 

is the k* row of 2° x 2" matrix B. Also, 

= [do di d 2 d 12 d )2 J 
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where the superscript T denotes matrix transpose. We note that B is a 2“ x 2“ square 
matrix and B = B T , hence B is symmetric. In other words, it is an orthogonal matrix 
[Ref. 4], Therefore: 



BB T = 2"/ 



(3.4) 



where I is 2° x 2" identity matrix. Thus, 



B l = ±B T 

2 n 



(3.5) 



and the solution of (3.3) is 



d = —B T x 

2 ” 



(3.6) 



A computer program listing is given in Appendix B to calculate the coefficients in 
d. The program requires the number of inputs k, the constraint length v, and the vector 
x. The elements of x are from the set {x}, x=l, 2, ... M. Thus, x is the signal 
assignment of the corresponding trellis diagram. 

Once the relation between the input stream {bj} and {x}, x=l, 2, ... M is known, 
the encoder is realized according to logic variables {aj. In other words, a natural 
TCM/MFSK mapper may be implemented if the analytic description of the encoder is 
available. In an orthogonal MFSK signal constellation, all signals have equal euclidean 
distance from each other as stated in (2.3). Thus, the distance is of no importance when 
the signals in x are assigned. Instead, care should be given to avoiding a catastrophic 
code (a limited number of errors causes an infinite number of errors at the Viterbi 
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decoder output) and intersymbol inference should be kept in mind. Another important 
concern may be the complexity of the encoder. If the three rules below are followed, 
generally one ends up with a good code: 

1 . All elements of x should be equiprobable. 

2. The first 2 m elements must consist of the elements of x=l, 2, M. 

3. The frequency differences must be kept constant giving the maximum signal 
difference to b,. 

These rules are derived from a computer search of many signal assignments. Rule 1 
ensures to keep the optimality of decoding which is maximum-likelihood and optimal for 
the equally likely messages. All three rules are also intended to avoid catastrophic 
conditions and reduce Doppler interference. If we end up with a consistent mapping rule, 
this ensures that TCM/MFSK is not catastrophic. Rule 3, with the other two rules, 
guarantees that adjacent transitions have the maximum frequency separation which 
reduces Doppler interference. 

A frequency difference corresponds to the absolute difference between any two 
elements of {x}, x=l, 2, ... M whose associated trellis transitions differ only one bit. 
For example, the frequency difference due to bj is 

8. = | x(b v b 2 , b. b k J -x(b x b 2 -b f ,b k J\ (3.7) 

When all 5j’s j = 1 ,2, . . .n are kept constant, d has a minimum number of nonzero elements 
and transmitter complexity is minimal [Ref. 6]. 
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A list of the codes for TCM/4-FSK and TCM/8-FSK with several constraint lengths 
is given at the end of this chapter. In signal assignment for 4-ary codes, the frequency 



differences are kept constant in a periodic fashion, as follows: 



5 , =2 , 8 2 = 1 , 8 3 - 2 , 5 4 = 1 



(3.8) 



On the other hand, for 8-ary codes, the following frequency differences are used. 



For TCM/8-FSK, k is 2 and the codes with constraint length v up to 8 are considered. 
Hence, the maximum n = 2 + 8 = 10. Because the channel waveform depends on 
n = k + v bits, TCM/8-FSK waveforms with v = 8 will be a function of 10 bits. If the 
constraint length v is taken to be 2, TCM/8-FSK waveforms are determined by 
n = k+ v = 4 bits. Thus, only 8 U 8 2 , 5 3 and 5 4 will be considered. 

Once we have an analytic expression for the code we may find the natural mapping 
rule in terms of the logic variables a,, a, G {0,1}. We note that a, = l /i( l-bj and that 
multiplication with b ( variables corresponds to an exclusive or function in logic variables 
a i . The natural mapping is the process of mapping the output of the convolutional encoder 
to one of the M signals according to a given rule. For example, Gray encoding is a kind 
of mapping with the rule that adjacent signals have a change of only one bit. The natural 
mapping is defined in logical terms, not in physical quantities; hence, the constants in 
analytic expression have no meaning in the natural mapping. 

The codes in logic variables corresponding to analytic expressions in Table 3.3 are 
listed in Table 3.4. The variables {zj}, j = 1, 2, ... m are the inputs to the natural 




(3.9) 
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frequency mapper. (Obviously, Zj’s can also be viewed as the outputs of the convolutional 
encoder.) The mapping rules for TCM/4-FSK and TCM/8-FSK are found by examining 
the signal assignment of each scheme and are given in Table 3.1 and Table 3.2, 
respectively. 



TABLE 3.1 TCM/4-FSK NATURAL MAPPING RULE 



Z, 


Z 2 


Frequency Slot 


1 o 


0 


1 


0 


1 


2 


1 


0 


3 


1 


1 


4 



The analytic expression and consequently the natural mapping rule is a function of 
signal assignment. Thus it is possible to find other codes by changing the signal 
assignment within the rules given previously. In other words, the code for a given 
constraint length, code rate, and modulation is not unique. 
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TABLE 3.2 TCM/8-FSK NATURAL MAPPING RULE 





Z 2 


Z3 


Frequency Slot 


0 


0 


0 


1 


0 


0 


1 


2 


0 


• 


0 


4 


0 


t 


1 


4 


• 


0 


0 


5 


1 


0 


1 


6 


1 


1 


0 


7 


1 


1 


1 


8 



TABLE 3.3 ANALYTIC DESCRIPTIONS OF THE CODES 



Code 


TCM/MFSK 


Constraint 


Analytic 


Rate (r) 




Length (v) 


Description 


1/2 


4-ary 


2 


x=2.5-0.5b 2 -b 1 b 3 


1/2 


4-ary 


3 


x=2.5-0.5b2b4-b,b 3 


1/2 


4-ary 


4 


x = 2 . 5-0 . 5b 2 b 4 -b 1 b 3 b 5 


1/2 


4-ary 


5 


x = 2 . 5 -0 . 5b 2 b 4 b 6 -b,b 3 b 5 


1/2 


4-ary 


6 


x=2.5-0.5b 2 b 4 b 6 -b 1 b 3 b 5 b 7 


1/2 


4-ary 


7 


x = 2 . 5-0 . 5b 2 b 4 b 6 b g -b,b 3 b 5 b 7 


2/3 


8-ary 


2 


x=4.5-0.5b 3 -b 2 -2b,b 4 


2/3 


8-ary 


3 


x=4.5-0.5b 3 -b 2 b J -2b 1 b 4 


2/3 


8-ary 


4 


x = 4 . 5 -0 . 5b 3 b 6 -b 2 bj-2b! b 4 


2/3 


4-ary 


5 


x = 4 . 5-0 . 5 b 3 b 6 -b 2 b 5 -2b 1 b 4 b 7 


2/3 


8-ary 


6 


x=4.5-0.5b 3 b6-b 2 bjb 8 -2b,b 4 b7 


2/3 


8-ary 


7 


x = 4. 5-0. 5b 3 b 6 b9-b 2 b5b8-2b,b 4 07 


2/3 


8-ary 


8 


x = 4 . 5 -0. 5030609 - 020503 - 20 , 04070,0 
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TABLE 3.4 CODES IN LOGIC VARIABLES 



Code Rate 
(r) 


TCM/MFSK 


Constraint 

Length 


Z, 


Z 2 


z. 


1/2 


4-ary 


2 


a!©a3 


a2 


- 


1/2 


4-ary 


3 


aj©a3 


a 2 ©a 4 


- 


1/2 


4-ary 


4 


a 1 ffia 3 ©a 5 


a 2 ©a 4 


- 


1/2 


4-ary 


5 


a 1 ©a 3 ©a 5 


a 2 ©a 4 ©a 6 


- 


1/2 


4-ary 


6 


a 1 ffia 3 ©a 5 ©a 

7 


a2 a 4 a 6 


- 


1/2 


4-ary 


7 


a 1 ©a 3 ©a 5 ©a 

7 


a 2 ©a 4 ©a 6 ©a 

8 


- 


2/3 


8 -ary 


2 


a j © a 4 


*2 


*3 


2/3 


8 -ary 


3 


a,©a 4 


a 2 ©a 5 


*3 


2/3 


8 -ary 


4 


a!©a 4 


a 2®a 5 




2/3 


8 -ary 




aj ©a 4 ©a 7 


a^a. 




2/3 


8 -ary 


6 


a, a 4 a 7 


a 2 ®a 3 ©a s 


a3©a5 


2/3 


8 -ary 


4 


a, a 4 ffi a^ 


a 2 ©a 5 ©a s 


a 3 ©a 6 ffia 9 J 


2/3 


8 -ary 


8 


a 1 ©a 4 ©a 7 ©a 

10 


a 2 ©a 3 ©a 8 


a 3 ©a ts ©a 9 1 



® denotes exclusive or 
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IV. ENCODER DESIGN EXAMPLES 



A. ENCODER DESIGN 1 

In this example, a rate Vi TCM/4-FSK encoder is considered. The constraint length 
v is chosen to be 3. There are v = 3 memory elements and T = 8 trellis states, each 
having 2 k = 2 transitions. Normalized squared free Euclidean distance is found by using 
the program in Appendix A, which is 4. The input to this program is the matrix N. The 
trellis structure, signal assignment, matrix N, and input/state variables are shown in 
Figure 4.1. 




1 2 
3 4 

5 6 



N= 



7 8 

1 2 



3 4 



5 6 



7 8 



Matrix N 




Trellis / Slgnsl Assignment Input / State Variables 
(The elements of x) 



Figure 4. 1 
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The rows of N are the labels of the destination states, respectively. For example, 
the first state has transitions to state 1 and state 2; hence, the first row of N has the 
labels 1 and 2. The second state has transitions to state 3 and state 4, and the elements 
of the second row are 3 and 4 and so on. 

The signal assignment is made by keeping the frequency differences constant; 
hence, 

6 , = 2 , 5 2 = 1 , 6 3 = 2 , 5 4 = 1 

For example, the transition from state 1 to state 1 occurs when the first frequency is 
transmitted. If the third frequency is transmitted then the new encoder state is the second 
one. The definition of frequency difference is stated in (3.7). The adjacent transitions are 
different in the b, bit only, thus the frequency difference 6, due to b, is 2. The signal 
assignment in Figure 4.1 is put in vector form as: 

x = [1 32431422413423 1] 

The input of the program given in Appendix B is the x vector. The program also requires 
the number of inputs k and the constraint length v. The outputs of this program are the 
coefficients in d which are found to be: 
d 0 = 2.5 
d i3 = -1 
d 24 — -0.5 

Thus, from (3.3) the analytic description of the code is 
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x = 2.5 -b l b 3 -0.5b 2 b i 




0.5 



Figure 4.2 Analytic Description Transmitter for TCM/4-FSK, 8 States 

The analytic description of the transmitter may be implemented as in Figure 4.2. The 
natural mapper implementation in terms of the logic variables a, is as follows. If we make 
the change of variables 

jc = 2.5 -y 1 -0.5y 2 

where y 2 = b,b 3 and y 2 = bp 4 , then since the multiplication in bj corresponds to an 
exclusive-or function in logic variables, we can find the new variables as 

z, = a x ®a 3 

z 2 = a 2 ®a A 
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These new variables z, t z 2 are input to the natural mapper. The mapping rule is found 
by examining the signal assignment. Thus, the natural mapper implementation may be 
as in Figure 4.3 with the following mapping rule. 



TABLE 4.1 



Zl 


z 2 


Frequency 


0 


0 


1 


0 


1 


2 


1 


0 


3 


1 


1 


4 




Zi 



MFSK 

Mapper 



sjt) 



fc 



Figure 4.3 Natural Mapper for r = x h, v = 3, TCM/4-FSK 
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The bandwidth efficiency of rate Vi TCM/4-FSK scheme can be found front (2.7) 



as 



k_ 

2 m 



- blsIHz 

4 



Uncoded 16-FSK has the same bandwidth efficiency. The ACG is 0 dB with respect to 
uncoded 16-FSK (Table 2.1). 



B. ENCODER DESIGN 2 

A rate 2/3 TCM/8-FSK with constraint length v = 4 is introduced in this example. 
There are 4 memory bits and 2 V = 16 trellis states. The number of transitions at each 
state is 2 k = 4. The squared free euclidean distance is found to be 3. The trellis diagram, 
signal assignment, matrix N, and input/state variables are shown in Figure 4.4. 

The bandwidth efficiency of the rate 2/3 TCM/8-FSK encoder is (2.7) 

— = — blsIHz 
2 m 4 

Uncoded 16-FSK has the same bandwidth efficiency. The ACG with respect to uncoded 
16-FSK is 1.76 dB (Table 2.1). 
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1 , 5 , 3 , 7 



2 , 4 , 4, 9 



5 , 1 , 7,3 



6 , 2 , 8 , 4 



3 , 7 , 1,5 



4 , 8 , 2, 4 



7 , 3 , 5,1 



8 , 4 , 4 , 2 



2 , 4 , 4 ,$ 



1 , 5 , 3 , 7 



4 , 2 , 8 , 4 



5 , 1 , 7 , 3 



4 . 8 . 2. 4 

3 . 7 . 1.5 



$, 4 , 4,2 



7 , 3 , 5,1 



N/ L 



l 






M 



\\m 



1 


2 


3 


4 1 


5 


4 


7 


8 


9 


10 


11 


12 


13 


14 


15 


14 


1 


2 


3 


4 


5 


4 


7 


8 


9 


10 


11 


12 


13 


14 


13 


14 


1 


2 


3 


4 


5 


4 


7 


8 


9 


10 


11 


12 


13 


14 


15 


14 


1 


2 


3 


4 


5 


4 


7 


8 


9 


10 


11 


12 


13 


14 


15 


14 



Trellis Metrlx 




Input /St ste VAxlmbles 



Trellis /Signal Assignment 



Figure 4.4 

When the signal assignment in Figure 4.4 is put in a vector x as 



x = [1537264851736284 
3715482673518462 
2648153762845173 
4826371584627351] 



and input to the program in Appendix B, the d coefficients are found to be 
do = 4.5 
d 14 = -2 
d 25 = -1 
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Hence, the analytical description of the set {x}, x = 1,2, 



8 in terms of b; is 



jc = 4.5 ~2b l b 4 -bp s -0.5 bp 6 



The above equation may be written as 

x = 4.5 -2y, -y 2 -0.5y 3 



where 

yi = bp 4 
y z = bp 5 
y 3 = bp 6 

Therefore, z, variables in terms of a ( are 

Zj = aj(Ba 4 

z 2 = a 2 ®a 5 
z 3 = a 3 ®a 6 

with the following mapping rule (Table 4.2). The natural mapper implementation of rate 
2/3 TCM/8-FSK with 16 states is shown in Figure 4.5. 
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TABLE 4.2 



2b 


*2 


Z ! 


Frequency 


0 


0 


0 


4 


• 


Q 


0 


2 


0 


1 


!) 


3 


t 


1 


0 


4 


0 


0 


1 


5 


1 


0 


1 


6 


0 


1 


1 


7 


1 


1 


1 


8 




fa 



Figure 4.5 Natural Frequency Mapper for r = 2/3, v = 4, TCM/8-FSK 
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C. ENCODER DESIGN 3 



A rate Vi TCM/4-FSK encoder with 32 states is considered. The constraint length 
is v = log 2 32 = 5. Each state has 2 transitions. The squared free Euclidean distance is 
6; ACG with respect to uncoded 16-FSK is 1.7 dB (Table 2.1). The matrix N which is 
the definition of the trellis diagram is given in (4.1) (following page). 

The signal assignment keeps the following frequency differences constant 

5i=2 , 5 2 = 1, 5 3 = 2, S 4 = l, 6 5 =2, S 6 = l 

The vector x is 

[13241422413423131421324 
jc= 423124132413423113243142 
4231241331421324] 

The above signal assignment yields the following analytic expression 

x = 2.5 -b 2 bjb 6 -b x b 3 b 5 

and the z, variables are 
Z\ ci j d j CD d$ 

Z 2 = d 2 ®d 4 ®d 6 

with the mapping rule given in Table 3.1. The natural mapper implementation is shown 
in Figure 4.6. 
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( 4 . 1 ) 
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Figure 4.6 Natural Frequency Mapper rate = Vi, v = 5, TCM/4-FSK 



D. ENCODER DESIGN 4 

A TCM/8-FSK encoder with 32 states is given to make the subject clearer. The 
rate is r = 2/3; hence k = 2, m = 3. The constraint length v is 5 and ACG is 1.76 dB. 
The matrix N is given in (4.2). 
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12 3 4 
5 6 7 8 
9 10 11 12 
13 14 15 16 
17 18 19 20 
21 22 23 24 
25 26 27 28 
29 30 31 32 
12 3 4 
5 6 7 8 
9 10 11 12 
13 14 15 16 
17 18 19 20 
21 22 23 24 
25 26 27 28 
29 30 31 32 
12 3 4 
5 6 7 8 
9 10 11 12 
13 14 15 16 
17 18 19 20 
21 22 23 24 
25 26 27 28 
29 30 31 32 



(4.2) 
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In this frequency assignment the following differences are kept constant (3.9) 

6j=4, S 2 = 2, 6 2 = 1, 5 4 =4, S 5 =2, S 6 = l, 8 7 = 4 

These constraints yield the following x vector. 

[15372648517362843715 
48267351846226481537 
62845173482637158462 
x = 73515173628415372648 

73518462371548266284 
51732648153784627351 
48263715 ] 

The x vector has 2 v+k = 128 elements. The following analytic expression is obtained by 
the signal assignment above 

x = 4.5 - 0.5 bp 6 - b 2 b 5 - 2bpp n 

This analytic expression corresponds the following z { variables with the rule given in 
Table 3.2 



h = a,®a 4 ®a 7 
z 2 = a 2 ®a 5 
z 3 = a 3 ®a 6 

Finally, the natural mapper can be implemented as shown in Figure 4.7. 
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Figure 4.7 Natural Frequency Mapper, r = 2/3, v = 5, TCM/8-FSK 
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V. DISCUSSION AND CONCLUSIONS 



The objective of this thesis was to investigate trellis coding of M-ary orthogonal 
signals. The power efficiency of this class of signaling is attractive to many 
communication systems when the environment is power-limited. An example of such an 
environment is the generation of low probability of intercept (LPI) signals. MFSK is used 
as an example of orthogonal signaling, but the discussion holds for other orthogonal 
signalings such as Pulse Position Modulation (PPM). 

A method of encoding orthogonal signals by using a smaller constellation with 
convolutional coding for a given bandwidth efficiency was considered. The free 
Euclidean distance of TCM/MFSK codes plays an important role when 
maximum-likelihood Viterbi soft-decision decoding is used on an AWGN channel. A 
computer program has been written to find the free Euclidean distance. A computer 
search of rate l /i TCM/4-FSK and rate 2/3 TCM/8-FSK shows that the normalized 
squared free Euclidean distance increases linearly with the constraint length, as; 

-*L = l-J + 1 

2 E L k J 

where v is the constraint length, k is the number of input bits, and [_*J * s the "floor" 
operation (the integer part). The normalized distance between any pair of signal points 
is taken to be unity; i.e., 
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d tj = J2E = 1 for all if ( i *j) 

The investigation shows that an arbitrarily small error rate can be achieved with a 
TCM/MFSK scheme by using a larger v. 

ACG comparisons are made with respect to uncoded 16-FSK having the same 
bandwidth efficiency. ACGs of 4.39 dB and 4.77 dB for TCM/4-FSK and TCM/8-FSK 
with constraint length 10 are observed, respectively. This means TCM/MFSK can 
maintain better performance while using the same bandwidth than uncoded MFSK. 

Convolutional coding with rates k/m are proposed to be used as the underlying 
code. This choice brings several advantages. First, efficient soft-decision Viterbi 
decoders can be used. Second, since rate k/m convolutional codes are used, the signal 
space is expanded as much as 2 m /2 k in terms of convolutional encoder input-output 
relations. It is shown that more coding gain is achievable with this method than the 
conventional approach. Doubling the signal space provides the greatest coding gain while 
requiring a minimum increase in system complexity [Ref. 7]. Finally, an outer 
coding structure may be used for further improvements in system performance. 

A computer program was written to find the analytic description of the codes. Once 
the analytic description is known, the natural mapper implementation is straightforward. 
The basic rules of frequency assignment for minimal transmitter complexity were 
discussed. The key is to keep the frequency differences constant. It is possible to fmd 
different codes by changing the signal assignment only. Encoder design examples are 
worked out which are also significant as another way of convolutional code design. 
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Although ACG gives an idea about the performance on an AWGN channel, the 
actual coding gain of trellis coded orthogonal signaling on both AWGN channels and 
fading channels for both coherent and noncoherent demodulation is required to complete 
our understanding of the performance that can be expected from this type of system 
under realistic conditions. 
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APPENDIX A 



%%%%%%%%%%%%% PROGRAM DISTANCE %%%%%%%%%%%% 



%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

% 

% This program is written in Matlab to find the free distance of a 
% trellis. The transitions from each state are saved in a matrix N. 

% This matrix is saved in a filename. m file. 

% The first row of the matrix N consists of the #'s of the states which 
% can be reached from state 1 in 1 transition .And the second row is for 
% the second state... (The first state is l,not 0) .The program takes care 
% of the rest. The parellel transtions are not considered. The program asks 
% the name of the input file. The output of the program is the squarred 
% free euclidean distance of the trellis described by the matrix N 
% 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

clear; 

% Input initilaziation 

r»input (' Enter the input file name ','s'); 

N=eval (r) ; 

% s : # of the states 

% t : # of the transitions at each state 
[s, t] = size (N) ; 

% Initilazation of the state transition matrix 
B=zeros (s, s) ; 
for i«l:s 

for j-l:t 

z-N(i, j) ; 

B (i, z) -z; 

end 

end 

for i*l:s 

for j-l:s 

if j*»B(i,j) 

B (i, j) =1; 
else 

B(i, j)-100; 
end; 

end; 

end; 

% Calculation of the shortest distances 
for k=l:s 

for i»l:s 

if B(i,k)~=100 
for j=l:s 

length=B ( i, k) +B (k, j ) ; 
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if length<B (i, j ) 

B (i, j) -length; 
end; 

end; 

end; 

end; 

end 

% Calculation of the free distance w.r.t. all zero codeword 
for j =1 : s 

if B(l,j)~ 1 

DF(1, j)-B(j,l)+B(l, j) ; 

DF (1 , 1) -100; 
else 

DF (1 , j ) -100; 

end 

end 

% Squarred free distance 
d2free-min (DF) 

% Pairwise Comparision for Checking 

% If the result of the pairwise comparision does not 

% agree the trellis setup is wrong 

dfree-100 ; 

yeni-B; 

for i-l;s 

if B (i, i) = = 1 
yeni (i, i) = 100; 
for j«l:s 

free-yeni (i, j) +yeni ( j , i) ; 
if free<dfree 
dfree=free; 

end 

end 

end 

end 

for i=l ; s 

for j=l:s 

if yeni (i, j) --1 
yeni (i , j ) -100; 
for x-l:s 

if yeni (i, x) --1 
for k-l:s 

if B ( j , k) ==B (x, k) 
free-B ( j , k) +1; 
if free<dfree 
dfree-free; 

end 

end 

end 

end 

end 

end 

end 

end 

dfree 
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APPENDIX B 



%%%%%%%% PROGRAM CODE %%%%%%%%% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

% This program finds the coefficients for the analytic desciption of a 
% TCM/MFSK scheme. The inputs to the program are the # of the inputs k, 

% the constraint length v, and a file. The input file contains a vector X 
% which is the signal assignment of the trellis. The number of the 
% elements is equal to the total number of the transitions. 

% 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

clear 

% Input initilazion 

m=input ( ' Enter the number of input bits # /'s'); 

r»input ( # Enter the costraint length ' 

t»input (' Enter the name of the input file ','s'); 

% k : # of the input bits 
% v : the constraint length 
% X : the signal vector 

% n : # of input bits plus the number of the memory bits 

k=eval (m) ; 
v=eval (r) ; 

X=eval ( t) ; 
n=k+v; 

% Creating Hadamard matrix 
Bshadamard (n) ; 

% Calculating coefficients 
D=B'*X'/2^n; 



% Determining the place of the coeffients 
d=f ind (D' ) 
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